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PRELIMINARY AMENDMENT 

Box Patent Application 

Assistant Commissioner for Patents 

Washington, D.C. 20231 

Dear Sir: 

Please enter the following preliminary amendment into the above-referenced 
application. 

ABSTRACT 

Please insert the attached abstract into the application as the last page 
thereof. 

SPECIFICATION 

Please enter the substitute specification submitted herewith as the 
specification to be examined. A clean copy of the substitute specification is included 
in Appendix B. A marked up copy of the substitute specification showing the 
amendments is included in Appendix A. 



CLAIMS 

Please delete claims 1-10 as follows. Please enter new claims 11-18 as 
follows. 

1 1 . (NEW) A method of channel estimation in a cellular communication 
system in which symbols are transmitted between a mobile station and a base 
station, the method comprising: 

receiving via a communication channel an incoming signal that incorporates symbols 
spread by a spreading code in a CDMA system, the incoming signal including data 
symbols and pilot symbols, a unique spreading code having been used respectively 
for pilot symbols and data symbols; 

providing a first estimate for the channel impulse response of that 
communication channel and utilising that first estimate to process the incoming 
signal to generate an estimated data symbol; 

using said estimated data symbol to generate a soft output feedback decision; 

and 

using the soft output feedback decision to make a further estimate of the 
channel impulse response. 

12. (NEW) A method according to claim 1 1 , wherein the channel impulse 
response is estimated by performing a correlation process between a conjugate 
version of the spreading code and the incoming signal, and multiplying the results of 
that process by the soft feedback decision value as a weighting value for the channel 
impulse response. 

13. (NEW) A method according to claim 1 1 , in which a long spreading code 
has additionally been used to spread the symbols prior to transmission. 



14. (NbW) r A method of channel estimation in a cellular communication 
system in which symbols are transmitted between a mobile station and a base 
station, the method comprising: 

receiving via a communication channel an incoming signal which incorporates 
symbols spread by a spreading code; 

providing a first estimate for the channel impulse response of that 
communication channel and utilising that first estimate to process the incoming 
signal to generate an estimated data symbol; 

making a further estimate of the channel impulse response by performing a 
correlation process between a conjugate version of the spreading code and the 
incoming signal, and multiplying the results of that process by a feedback decision 
based on said estimated data symbol as a weighting value for the channel impulse 
response. 

1 5. (NEW) A method according to claim 14, wherein the feedback decision 
is a soft output. 

16. (NEW) A method according to claim 14, wherein the feedback decision 
is a hard output. 

17. (NEW) Circuitry for estimating a channel impulse response in a cellular 
communication system in which symbols are transmitted between a mobile station 
and a base station, the circuitry comprising; 

a receiver for receiving via a communication channel an in coming signal 
which incorporates symbols spread by a spreading code; 

a channel impulse response generator for providing a first estimate for the 
channel impulse response of that communication channel; 



a data symbol generator for generating an estimated data symbol using the 
first estimate of the channel impulse response; 

circuitry for providing a further estimate of the channel impulse response by 
performing a correlation process between a conjugate version of the spreading code 
and the incoming signal, and multiplying the results of that process by a feedback 
decision based on said estimated data symbol as a weighting value for the further 
estimate of the channel impulse response. 

1 8. (NEW) Circuitry for estimating a channel in a cellular communication 
system in which symbols are transmitted between a mobile station and a base 
station, the circuitry comprising: 

a receiver for receiving an incoming signal via a communication channel, the 
incoming signal incorporating symbols spread by a spreading code in a CDMA 
system, the incoming signal including data symbols and pilot symbols, a unique 
spreading code having been used respectively for pilot symbols and data symbols; 

a channel impulse response generator for providing a first estimate for the 
channel impulse response of that communication channel; 

a data symbols generator for generating an estimated data symbol using that 
first estimate; 

a soft output feedback decision generator for using said estimated data 
symbol to generate a soft output feedback decision; 

wherein the soft output feedback decision is used to make a further estimate 
of the channel impulse response. 



REMARKS 

The above preliminary amendment is made to insert an abstract page into the 
application, to enter a substitute specification, to delete claims 1-10 and to add new 
claims 1 1-18 for examination. 

Applicant respectfully requests that this preliminary amendment be entered 
into the record prior to calculation of the filing fee and prior to examination and 
consideration of the above-identified application. 

If a telephone conference would be helpful in resolving any issues concerning 
this communication, please contact Applicant's attorney of record, Michael B. Lasky 
at 952-912-0527. 
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Appendix A 

Marked Up Version of the Sustitute Specification 

The present invention relates to a method and circuitry for channel estimation in a 
cellular communication system. 

The invention is particularly concerned with a cellular communication system in 
which data is transmitted as a plurality of data symbols over a sequence of time 
slots. As is known, in a CDMA system, data is encoded for transmission by 
modulating data symbols to be transmitted by a unique spreading code for each of a 
plurality of communication channels. Within each cell of a cellular communication 
system, spreading codes allow for a plurality of different mobile stations to 
communicate with a base station on selectively coded channels. 

When a signal is transmitted between a base station and a mobile station (either on 
the uplink or the downlink), the signal receiving unit needs to establish from the 
signal which it has received some information about the communication path along 
which the signal has travelled. This is referred to herein as "channel estimation" and 
is carried out in a channel estimation unit which generates a channel impulse 
response. Various techniques are known for channel estimation. The channel 
impulse response is required in order to properly decode and demodulate incoming 
data. 

In an article in the IEEE Journal a Selected Areas in Communications entitled "A 
Soft-Output Bidirectional Decision Feedback Equilization Technique for TDMA" (II 
(1993) September, No 7) a bidirectional adaptive equilization technique that takes 



advantage of the periodic transmission of training bits in the TDMA frame structure to 
estimate the location of deep fade within a time slot is proposed. 

US-A-5,732,114 discloses a method of detecting reference signals for a digital data 
receiver utilising a correlation process between a conjugate of a reference symbol 
and an incoming signal. 

In an earlier proposed CDMA system explained in the Technical Report of IEICE 
(1996 - 02) entitled "Variable Rate Data Transmission on Single Code Channels in 
DS/CDMA" by Okomura et a/, data is transmitted at variable transmission rates. 
Channel estimation is carried out from pilot symbols spaced periodically in the data 
stream. 

One difficulty associated with this proposal is that only the energy from the pilot 
symbols is available for channel estimation, and the signal to noise ratio of pilot 
symbols tends to be poor. Thus, filtering is required to get more reliable estimates. 
In particular, in the case of mobile stations moving at high speeds, it may be difficult 
to apply a suitable filtering technique due to short channel coherence time. 

It is therefore advantageous to use data symbols in addition to the pilot symbols for 
channel estimation. The present invention seeks to provide a technique for 
improving channel estimation using estimated data symbols. 

According to one aspect of the invention there is provided a 



method of channel estimation in a cellular communication system in which symbols 
are transmitted between a mobile station and a base station, the method comprising: 
receiving via a communication channel an incoming signal that incorporates symbols 
spread bv a spreading code in a CDMA system, the i ncoming signal including data 
symbols and pilot symbols, a unique spreading code h aving been used respectively 
for pilot symbols and data symbols: 

providing a first estimate for the channel impulse response of that 
communication channel and utilising that first estimate to process the incoming 
signal to generate an estimated data symbol; 

using said estimated data symbol to generate a soft output feedback decision; 

and 

using the soft output feedback decision to make a further estimate of the 
channel impulse response. 

[The invention is particularly but not exclusively applicable to a CDMA system in 
which the incoming signal incorporates symbols spread by a spreading code.] 

According to another aspect of the invention there is provided a method of channel 
estimation in a cellular communication system in which symbols are transmitted 
between a mobile station and a base station, the method comprising: 

receiving via a communication channel an incoming signal which incorporates 
symbols spread by a spreading code; 



providing a first estimate for the channel impulse response of that 
communication channel and utilising that first estimate to process the incoming 
signal to generate an estimated data symbol; 

making a further estimate of the channel impulse response by performing a 
correlation process between a conjugate version of the spreading code and the 
incoming signal, and multiplying the results of that process by a feedback decision 
based on said estimated data symbol as a weighting value for the further estimate of 
the channel impulse response. 

The incoming signal can include data symbols and pilot symbols, with a unique 
spreading code having been used respectively for pilot symbols and data symbols. 

According to the further aspect of the invention, the feedback decision may be a soft 
output or a hard output. 

A further aspect of the invention provides circuitry for estimating a channel impulse 
response in a cellular communication system in which symbols are transmitted 
between a mobile station and a base station, the circuitry comprising: 

a receiver for receiving via a communication channel an in coming signal 
which incorporates symbols spread by a spreading code; 

a channel impulse response generator for providing a first estimate for the 
channel impulse response of that communication channel; 

a data symbol generator for generating an estimated data symbol using the 
first estimate of the channel impulse response; 



circuitry for providing a further estimate of the channel impulse response by 
performing a correlation process between a conjugate version of the spreading code 
and the incoming signal, and multiplying the results of that process by a feedback 
decision based on said estimated data symbol as a weighting value for the further 
estimate of the channel impulse response. 

A still further aspect of the present invention provides circuitry for estimating a 
channel in a cellular communication system in which symbols are transmitted 
between a mobile station and a base station, the circuitry comprising: 

a receiver for receiving an incoming signal via a communication channel, the 
incoming signal incorporating symbols spread bv a spreading code in a CDMA 
system, the incoming signal including data symbols and pilot symbols, a unique 
spreading code having been used respectively for pi lot symbols and data symbols; 

a channel impulse response generator for providing a first estimate for the 
channel impulse response of that communication channel; 

a data symbols generator for generating an estimated data symbol using that 
first estimate; 

a soft output feedback decision generator for using said estimated data 
symbol to generate a soft output feedback decision; 

wherein the soft output feedback decision is used to make a further estimate 
of the channel impulse response. 

The invention also contemplates a base station and/or a mobile station incorporating 
circuitry according to the first aspects of the invention. 



For a better understanding of the present invention and to show how the same may 
be carried into effect, reference will now be made by way of example to the 
accompanying drawings in which: 

Figure 1 is a block diagram illustrating schematically symbol estimation in a 
basic CDMA system; 

Figure 2 is a block diagram illustrating a cellular communication system; 
Figure 3 is a diagram illustrating a transmission time slot structure; 
Figure 4 is a block diagram illustrating channel estimation with feedback; 
Figure 5 is a block diagram of circuitry in a channel impulse response unit; 
Figure 6 is a vector space diagram; and 

Figures 7 and 8 are simulations of acquisition probability against signal 
energy per bit to noise ratio. 

Figure 1 is a block diagram of a basic receiving system in a CDMA cellular 
communication network. Figure 2 is a block diagram 



illustrating a context in which the present invention may be used. That is, a CDMA 
mobile communication system allows a plurality of mobile stations MS1, MS2, MS3 
to communication with a base station BTS in a common cell via RF signals on 
respective channels CH1, CH2, CH3, These channels are distinguished from one 
another by the use of spreading codes in a manner which is known per se. 

Reverting now to Figure 1 , signals incoming at an antenna 30 are received by an RF 
unit 28 and supplied to an analogue-to-digital (A/D) converter 32. The A/D converter 
takes samples of the incoming signals at intervals i to generate digital sampled 
values x. It will readily be understood that a signal may arrive at the receiving 
circuitry having experienced multi-paths with differing propagation delays. Each 
digital sampled value Xj may include therefore components of a number of different 
symbols depending on the multipath effect. 

The digital signal is supplied from the A/D converter 32 to a channel impulse 
response unit 37 which provides an estimated channel impulse response H for use in 
processing the incoming signal. A reference code generator 22 supplies a reference 
code for despreading the incoming signal both to the channel impulse response unit 
37 and to a rake receiver 36. The reference code is denoted C and is the unique 
spreading code for the communication channel on which the signal is received. 
Thus, it matches the code which was used to spread the signal prior to transmission 
at the base station or mobile station respectively. Although Figure 1 illustrates 
circuitry located at each mobile station, similar circuitry may also be utilised at a base 
station. 



The rake receiver 36 additionally receives the input digital signal x and the estimated 
channel impulse response H. It generates an estimated transmitted signal u which is 
supplied to a decoder 31 . The decoder estimates symbol values s from the 
estimated signal u which has been derived from the actual incoming signal x, in a 
manner which is know per se. 

Operation of the rake receiver is known to a person skilled in the art but in any event 
is described in more detail in Application W096/24988 in the name of Nokia Mobile 
Phones Limited. 

The form in which data is transmitted is illustrated by way of example in Figure 3, 
The data symbols are transmitted in a sequence of time slots, each time slot 
including pilot symbols (PS) and coded data. Data is transmitted as a set of 
analogue symbols which are produced from originating digital data by i,q 
modulation. In a CDMA system the symbols are spread by a spreading code prior to 
transmission. In a typical system, each time slot has a duration of .625ms. The pilot 
symbols are introduced at the beginning and end of each time slot in a coherent 
system. These symbols are readily recognisable and so the beginning and end of 
each time slot can be identified for synchronisation purposes. In the circuit of Figure 
1, the pilot symbols are additionally identified and used in the estimation of the 
channel impulse response by the channel impulse response 37. 



Figure 4 illustrates an improved version of the circuitry of Figure 1 . In Figure 4, the 
estimated symbol s is fed back to the channel impulse response unit 37 via a 
feedback decision generator 20 which generates a feedback decision D from the 
estimated symbol s. That feedback decision D is used to improve the estimate of 
channel impulse response generated by the CIR unit 37. An SIR block 29 receives 
the incoming signal samples x and the estimated signal u and calculates the signal 
to information ratio used to generate noise variance a 2 for the feedback decision 
generator 20. 

Figure 5 illustrates how the feedback decision D is used to improve the estimated 
channel impulse response. Figure 5 illustrates a CIR unit 37 which is capable of 
providing an estimated channel impulse response H using both fedback estimated 
data symbols and incoming pilot symbols. The circuit comprises first and second 
matched filters 10,12 which are each used to apply the relevant despreading code to 
the incoming digital signal Xj n . The reference code generator 22 generates three 
codes. A spreading long code LC is applied to both of the match filters 10,12. A 
spreading short code W p for pilot symbols is applied only to the first match filter 10 
and a spreading short code W d for data is applied only to the second match filter 12. 
The output signals 14,16 from the first and second match filters 10,12 respectively 
are applied to an averaging circuit 18. The averaging circuit receives the decision 
feedback D and applies that to the signals 14,16 from the matched filters to generate 
the estimated impulse response H. 



The averaging circuit 18 includes a memory 40 for holding previous channel impulse 
responses and a weighting/averaging unit 42 for taking an average of the earlier 
channel impulse responses and applying a weighting function according to the 
fedback decision value D. 

The correlation process which is used to establish the estimate of channel impulse 
response can be expressed as: 

H n (t)= I D n *C* (i)x(i+t) =D n £ C*(i)x(i+t) (Equation 1) 

i i 

In the above equation, i is the sampling time index along the received signal, t is the 
resulting sampling time index of the correlation process, x(i+t) is the received data 
sample for that index, C* is the conjugate of the spreading code and D n is the n m 
decision feedback where n is the time index of the impulse response. It will be 
apparent that the decision feedback D may be in respect of pilot symbols and/or data 
symbols. 

In use, an initial estimate for the channel impulse response H would be established 
by the CIR unit 37 on the basis only of known incoming pilot symbols in a particular 
time slot in a manner which is known perse. Alternatively, an initial estimate can be 
made using Equation (1) without the weighting function D. Once a preliminary 
estimate of the channel impulse response has been established, this is held in the 
memory 40 and supplied to the rake receiver 36 for subsequent signal processing of 
data symbols in that time slot. The first data symbol to be estimated is then subject 



to feedback to prbvide a feedback decision D so that a new estimate for the channel 
impulse response H can be determined using the procedure outlined above and as 
illustrated in Figure 5. The spreading code conjugate C* in equation 1 is derived 
from merging the conjugates of Wd and LC (for data) or W p and LC (for pilot 
symbols). Then, an average is taken of the new estimate and the preliminary 
estimate. Thus, at each estimated symbol an improved channel impulse response is 
obtained for subsequent processing. The average can be taken coherently within 
any coherent period, that is part of a slot, one whole slot or a multiple of slots. 
Alternatively, a non-coherent average can be taken over any period. 

It will be apparent from Equation 1 that the correlation process can effectively be 
implemented without decision feedback as if all the transmitted data symbols are 
ones. The transmission feedback D acts as an average weighting factor to 
determine the estimated channel impulse response H. 

Another aspect of the implementation discussed herein is to determine the nature of 
the decision feedback D generated by the feedback decision generator 20. Figure 6 
illustrates the theory underlying one proposal which is to utilise a soft symbol 
decision feedback (SDFB) process by generating a soft feedback decision value D. 
Figure 6 illustrates a vector diagram of signal space. Vector A is a transmitted data 
signal. Vector A 1 is the received data. Vector A is the actual value we expect from 
the transmitted data signal since the rake receiver cannot be perfect. The received 
data symbol can be written as: 



A-a+n ' 
where n is the noise vector. 



(Equation 2) 



It is natural to assume that both the signal and the noise are zero mean. We have: 

E{a} =0 
E{n} =0 

E{A*}=E{a+n} =0 (Equation 3) 

In a binary (+1/-1) environment, we have also: 
E{AA*} =1 

E{A'A'*} =E{(a+n)(a+n)*} 
=E{a a*}+ E{nn*} 

=E{a a*}+ a 2 (Equation 4) 

where a 2 is the noise variance. 

With these statistics of the signal, the soft output can be adjusted to satisfy the 
conditions for feedback if the noise variance is known beforehand. For a simple 
solution, we suppose that the signal loss is minimum and can be neglected: 



U-a| -> o 



Then equation (3) becomes: 



E{A'A'*} =E{a a*}+ cr 



.2 



=E{A A*}+ a 2 



(Equation 5) 



Using this analysis, the feedback decision generator can generate a soft value D 
satisfying the following requirements: 



As an alternative, a hard decision feedback may be utilised by performing a binary 
cut-off, that is to provide either +1 or -1 to the input of the channel impulse response 
unit for each estimated symbol S. 

The effectiveness of both solutions is illustrated in Figures 7 and 8. In figures 7 and 
8, the x axis represents signal energy per bit to noise ratio (E b /N 0 ) and the y axis 
represents the scale of acquisition probability - that is the likelihood that all signal 
propagation paths have been captured in the received signal as estimated using the 
C.I.R. Figure 7 is the simulation results with the best possible channel estimation, 
and Figure 8 is the simulation results in which the channel information of a slot is 



E{D} = 0 



(Equation 6) 



E{DD*}= 1+a 2 



(Equation 7) 



obtained only frorn the pilot symbols in the same slot. These simulations have been 
generated based on the MATLAB program included as Annexe 1 . The simulation 
results compare the acquisition probabilities of the impulse response of different 
approaches of the decision feedback compared with no decision feedback. In both 
cases, it can be seen that the best results are given by a soft decision feedback for 
both pilot and data symbols. This is the top line in both simulations. The next line is 
for soft feedback decisions relating only to data symbols, that is without pilot symbol 
feedback being taken into account. 

The next line is for hard data feedback decisions in respect of data symbols only. 
The lowest line is where there is no decision feedback. 



ANNEXE 1 



MAT LAB program 

% Soft Decision Feedback 

% ARIB specification 

% QZD/NTC - FEB. 98 - Oulu 

clear, rand ( ' seed 1 , 0) , randn ( 1 seed' , 0) 

Datjtfeight = 0.5; SNR = -30:3:0; Fingers = 10; 

h=[0.7+j*0.1 / 0 / 0,0 / -0.5-j*0.3 / 0,0 / 0 / 0 / 0.3+j*0.5 / 0,0 / ... 

0.1-j*0.13,0,0,0,0,0,0,0.3-j*0.1,0,0,0,0.5+j*0.5, ... 

0, 0, 0, -0. 5-j ♦0.3,0,0,0,0,0. 3+ j*0. 5,0, 0,0. l-j*0. 13,... 

0,0,0,0, 0,0, 0.3-j*0.1] ; 
h=h/sqrt (h*h' ) ; 

AcqError = zeros (size (SNR, 2) , 4) ; 

for kkk = 1: size (SNR, 2) 

snr = SNR (kkk) % SNR in dB 

for nnn=l:1000 
DAT = (rand(l,40) >0.5) *2-l; 
% CODE 

CL_I = (rand (1,2560) >0.5) *2-l; CL_Q= (rand (1 , 2 560 ) >0.5) *2-l; 
C = CL_I + j*CL_Q; 

w=l;for k=l:6;w=[w w;w -w] ; end; w0=w (1 , : ) ; wl=w (2 , : ) ; clear W; 
W0=[];W1=[]; for k=l:40, W0= [W0 , wO] ; Wl= [Wl , wl] ; end 
CO = W0 .*C; CI = Wl . *C; 
% Spreading 

D=[] ; for k=l:40, D- [D , DAT (k) *wl] ; end % Data 
P=ones (size (D) ) ; % PILOT 

xl = (D + j*P) .* C; 
% Channel 

x2 = conv(xl,h) /2; 
% Noise 

nA=l/l0 x (snr/20) ; f or k=l : size (x2 , 2) 

n (k) =nA* randn* exp ( j *rand*2*pi) ;end 
x = x2 + n; 

%%%%%%%%%%%%%%%%%%% MF %%%%%%%%%%%%%%%%%%% 
IRPilot = zeros (39, 64); IRData = zeros (39, 64); 
for k=l:39, l=(k-l)*64; 
for m = 1 : 64 

IRPilot (k,m) = IRPilot (k,m) + x (l+m-l+ (1 : 64) ) *C0 (1+ (1 : 64) ) » ; 
IRData (k,m) = IRData (k,m) + x {l+m-l+ (1 : 64) ) *C1 (1+ (1 : 64) ) ' ; 
end 

end 

IRPilot = -j*IRPilot/64; 
IRPilotHard=sum (IRPilot (1:28, :) ) /28; 
IRData = IRData/64; 

%%%%%%%%%%%%%%%%%%% H' %%%%%%%%%%%%%%%%%%%%%%%%%% 
a = IRPilotHard .* con j (IRPilotHard) ; 
H = zeros (size (h) ) ,- 
for k = 1: Fingers, findf = 0; 
while ( findf < 1) 

[n,m] =max (a) ; a (m) =0 ; 



if (m<=si'ze (h, 2) ) , H (m) =IRPilotHard (m) ; findf=2; end 

end 

end 

%%%%%%%%%%%%%%%%%%% decision %%%%%%%%%%%%%%%%%%% 

hl=find(abs (h) >0.01) ; 

h2=hl; % Idea 

%h2=f ind (abs (H) >0 . 01) ; % blind 

PilotSoft = zeros(l,28); DataSoft = zeros(l,40); 
for n=l: Fingers 

a = x(h2 (n) ;size(C0,2) +h2 (n) -1) .* conj (CO) * conj (H(h2 (n) ) ) ; 
for k=l:28 

PilotSoft (k) = PilotSoft (k) + imag (sum(a ( (k-1) *64 + l :k*64) ) ) ; 
end 

a = x(h2(n) : size (CI, 2) +h2 (n) -1) .* conj (CI) * conj (H (h2 (n) ) ) ; 
for k=l:40 

DataSoft (k) = DataSoft (k) + real (sum(a ( (k-1) *64+l :k*64) )) ; 
end 

end 

PilotSoft = PilotSoft/64; 
DataSoft = DataSof t/64; 
DataHard = sign (DataSof t) ; 

5- 9- 9- 9- 9- 2- 9- 2- 2- 2- 9- 9- 2- 9- 9- 2- 9- 9- 5- AiTQ-rnrri 9- 9- 9- 2- 9- 9- 9* 9- 9- 9- 2- 9- 2- 9- 9- 9- 9- 2- 9* 

"o'o'o'o'o'o'o'o'^'o'o'o'o'b'o'o'3'o'o /iVCXcty-LIiy "5o*bo'ooo*5oooooo*6oooo 

a = 0; b = 0; C = 0; 

IRHard=zeros (1, 64) ; IRSof t=zeros (1 , 64) ; IRSof tp=zeros (1, 64) ; 
IRSof ter=zeros (1, 64) ; IRSof terp=zeros (1, 64) ; 
for k=l:39 

IRHard=IRHard+IRData(k, : ) *DataHard(k) ; a=a+abs (DataHard (k) ) ; 
IRSof t=IRSof t+IRData (k, : ) *DataSof t (k) ; b=b+abs (DataSoft (k) ) ; 

end 

for k=l:28 

IRSof tp= IRSof tp+IRPi lot (k, : ) *PilotSof t (k) ; c=c+abs 
(PilotSoft (k) ) ; 

end 

IRSoftp = IRSoft*Dat_Weight/b + IRSoftp * ( l-Dat_Weight) /c ; 

IRHard = IRHard*Dat_Weight/a + IRPilotHard* (l-Dat_Weight) ; 

IRSoft = IRSoft*Dat_Weight/b + IRPilotHard* (l-Dat_Weight) ; 

s> *********** * accj probability ************ 

a = IRPilotHard .* conj (IRPilotHard) ; 

for k=l : size (hi , 2 ) , [n, m] =max (a) ;h2 (k) =m; a (m) =0 ;end 

c=0; for k=l : size (hi, 2) , f or n=l : size (hi, 2) 

if (h2 (k) ==hl (n) ) , c=c+l; h2(k)=0; end 
end, end 
AErrNoDFB = size (hi, 2) -c; 
a = IRHard . * conj (IRHard) ; 

for k-l:size(hl,2) , [n,m] =max(a) ;h2 (k) =m;a(m) =0;end 
c=0; for k=l:size(hl,2) , f or n=l : size (hi , 2 ) 

if (h2 <k) ==hl (n) ) , c=c+l; h2(k)=0; end 
end , end 
AErrHardDFB = size (hi , 2 ) -c ; 
a = IRSoft .* conj (IRSoft) ; 

for k=l : size (hi, 2) , [n,m] =max(a) ;h2 (k) =m;a (m) =0;end 
c=0; for k=l : size (hi, 2) , for n=l : size (hi , 2) 



if (hi (k) ==hl (n) ) , c=c+l; h2(k)=0; end 
end, end 
AErrSof tDFB = size (hi, 2) -c; 
a = IRSoftp . * conj (IRSoftp) ; 

for k=l :size (hi, 2) , [n,m] =max(a) ;h2 (k) =m; a (m) =0 ; end 
c=0; for k=l : size (hi, 2) , for n=l : size (hi, 2) 

if (h2 (k) ==hl (n) ) , c=c+l; h2(k)=0; end 
end , end 
AErrSof tPDFB = size (hi, 2) -c; 
AcqError (kkk, : ) =AcqError (kkk, :)+... 
[AErrNoDFB , AErrHardDFB , AErrSof tDFB , AErrSof tPDFB] 
end % nnn 
end % kkk 



Appendix B 
Clean Version of the Substitute Specification 

The present invention relates to a method and circuitry for channel estimation in a 
cellular communication system. 

The invention is particularly concerned with a cellular communication system in 
which data is transmitted as a plurality of data symbols over a sequence of time 
slots. As is known, in a CDMA system, data is encoded for transmission by 
modulating data symbols to be transmitted by a unique spreading code for each of a 
plurality of communication channels. Within each cell of a cellular communication 
system, spreading codes allow for a plurality of different mobile stations to 
communicate with a base station on selectively coded channels. 

When a signal is transmitted between a base station and a mobile station (either on 
the uplink or the downlink), the signal receiving unit needs to establish from the 
signal which it has received some information about the communication path along 
which the signal has travelled. This is referred to herein as "channel estimation" and 
is carried out in a channel estimation unit which generates a channel impulse 
response. Various techniques are known for channel estimation. The channel 
impulse response is required in order to properly decode and demodulate incoming 
data. 

In an article in the IEEE Journal a Selected Areas in Communications entitled "A 
Soft-Output Bidirectional Decision Feedback Equilization Technique for TDMA" (II 
(1993) September, No 7) a bidirectional adaptive equilization technique that takes 



advantage of the periodic transmission of training bits in the TDMA frame structure to 
estimate the location of deep fade within a time slot is proposed. 

US-A-5,732,1 14 discloses a method of detecting reference signals for a digital data 
receiver utilising a correlation process between a conjugate of a reference symbol 
and an incoming signal. 

In an earlier proposed CDMA system explained in the Technical Report of IEICE 
(1996 - 02) entitled "Variable Rate Data Transmission on Single Code Channels in 
DS/CDMA" by Okomura et al, data is transmitted at variable transmission rates. 
Channel estimation is carried out from pilot symbols spaced periodically in the data 
stream. 

One difficulty associated with this proposal is that only the energy from the pilot 
symbols is available for channel estimation, and the signal to noise ratio of pilot 
symbols tends to be poor. Thus, filtering is required to get more reliable estimates. 
In particular, in the case of mobile stations moving at high speeds, it may be difficult 
to apply a suitable filtering technique due to short channel coherence time. 

It is therefore advantageous to use data symbols in addition to the pilot symbols for 
channel estimation. The present invention seeks to provide a technique for 
improving channel estimation using estimated data symbols. 

According to one aspect of the invention there is provided a 



method of channel estimation in a cellular communication system in which symbols 
are transmitted between a mobile station and a base station, the method comprising: 
receiving via a communication channel an incoming signal that incorporates symbols 
spread by a spreading code in a CDMA system, the incoming signal including data 
symbols and pilot symbols, a unique spreading code having been used respectively 
for pilot symbols and data symbols; 

providing a first estimate for the channel impulse response of that 
communication channel and utilising that first estimate to process the incoming 
signal to generate an estimated data symbol; 

using said estimated data symbol to generate a soft output feedback decision; 

and 

using the soft output feedback decision to make a further estimate of the 
channel impulse response. 

According to another aspect of the invention there is provided a method of channel 
estimation in a cellular communication system in which symbols are transmitted 
between a mobile station and a base station, the method comprising: 

receiving via a communication channel an incoming signal which incorporates 
symbols spread by a spreading code; 

providing a first estimate for the channel impulse response of that 
communication channel and utilising that first estimate to process the incoming 
signal to generate an estimated data symbol; 

making a further estimate of the channel impulse response by performing a 
correlation process between a conjugate version of the spreading code and the 



incoming signal, and multiplying the results of that process by a feedback decision 
based on said estimated data symbol as a weighting value for the further estimate of 
the channel impulse response. 

The incoming signal can include data symbols and pilot symbols, with a unique 
spreading code having been used respectively for pilot symbols and data symbols. 

According to the further aspect of the invention, the feedback decision may be a soft 

i 

output or a hard output. 

A further aspect of the invention provides circuitry for estimating a channel impulse 
response in a cellular communication system in which symbols are transmitted 
between a mobile station and a base station, the circuitry comprising: 

a receiver for receiving via a communication channel an in coming signal 
which incorporates symbols spread by a spreading code; 

a channel impulse response generator for providing a first estimate for the 
channel impulse response of that communication channel; 

a data symbol generator for generating an estimated data symbol using the 
first estimate of the channel impulse response; 

circuitry for providing a further estimate of the channel impulse response by 
performing a correlation process between a conjugate version of the spreading code 
and the incoming signal, and multiplying the results of that process by a feedback 
decision based on said estimated data symbol as a weighting value for the further 
estimate of the channel impulse response. 



A still further aspect of the present invention provides circuitry for estimating a 
channel in a cellular communication system in which symbols are transmitted 
between a mobile station and a base station, the circuitry comprising: 

a receiver for receiving an incoming signal via a communication channel, the 
incoming signal incorporating symbols spread by a spreading code in a CDMA 
system, the incoming signal including data symbols and pilot symbols, a unique 
spreading code having been used respectively for pilot symbols and data symbols; 

a channel impulse response generator for providing a first estimate for the 
channel impulse response of that communication channel; 

a data symbols generator for generating an estimated data symbol using that 
first estimate; 

a soft output feedback decision generator for using said estimated data 
symbol to generate a soft output feedback decision; 

wherein the soft output feedback decision is used to make a further estimate 
of the channel impulse response. 

The invention also contemplates a base station and/or a mobile station incorporating 
circuitry according to the first aspects of the invention. 

For a better understanding of the present invention and to show how the same may 
be carried into effect, reference will now be made by way of example to the 
accompanying drawings in which: 



Figure 1 is a block diagram illustrating schematically symbol estimation in a 

basic CDMA system; 

Figure 2 is a block diagram illustrating a cellular 
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communication system; 

Figure 3 is a diagram illustrating a transmission time slot structure; 
Figure 4 is a block diagram illustrating channel estimation with feedback; 
Figure 5 is a block diagram of circuitry in a channel impulse response unit; 
Figure 6 is a vector space diagram; and 

Figures 7 and 8 are simulations of acquisition probability against signal 
energy per bit to noise ratio. 

Figure 1 is a block diagram of a basic receiving system in a CDMA cellular 
communication network. Figure 2 is a block diagram 



illustrating a context in which the present invention may be used. That is, a CDMA 
mobile communication system allows a plurality of mobile stations MS1 , MS2, MS3 
to communication with a base station BTS in a common cell via RF signals on 
respective channels CH1 , CH2, CH3. These channels are distinguished from one 
another by the use of spreading codes in a manner which is known per se. 

Reverting now to Figure 1 , signals incoming at an antenna 30 are received by an RF 
unit 28 and supplied to an analogue-to-digital (A/D) converter 32. The A/D converter 
takes samples of the incoming signals at intervals i to generate digital sampled 
values x. It will readily be understood that a signal may arrive at the receiving 
circuitry having experienced multi-paths with differing propagation delays. Each 
digital sampled value Xj may include therefore components of a number of different 
symbols depending on the multipath effect. 

The digital signal is supplied from the A/D converter 32 to a channel impulse 
response unit 37 which provides an estimated channel impulse response H for use in 
processing the incoming signal. A reference code generator 22 supplies a reference 
code for despreading the incoming signal both to the channel impulse response unit 
37 and to a rake receiver 36. The reference code is denoted C and is the unique 
spreading code for the communication channel on which the signal is received. 
Thus, it matches the code which was used to spread the signal prior to transmission 
at the base station or mobile station respectively. Although Figure 1 illustrates 
circuitry located at each mobile station, similar circuitry may also be utilised at a base 
station. 



The rake receiver 36 additionally receives the input digital signal x and the estimated 
channel impulse response H. It generates an estimated transmitted signal u which is 
supplied to a decoder 31 . The decoder estimates symbol values s from the 
estimated signal u which has been derived from the actual incoming signal x, in a 
manner which is know perse. 

Operation of the rake receiver is known to a person skilled in the art but in any event 
is described in more detail in Application W096/24988 in the name of Nokia Mobile 
Phones Limited. 

The form in which data is transmitted is illustrated by way of example in Figure 3. 
The data symbols are transmitted in a sequence of time slots, each time slot 
including pilot symbols (PS) and coded data. Data is transmitted as a set of 
analogue symbols which are produced from originating digital data by i,q 
modulation. In a CDMA system the symbols are spread by a spreading code prior to 
transmission. In a typical system, each time slot has a duration of .625ms. The pilot 
symbols are introduced at the beginning and end of each time slot in a coherent 
system. These symbols are readily recognisable and so the beginning and end of 
each time slot can be identified for synchronisation purposes. In the circuit of Figure 
1 , the pilot symbols are additionally identified and used in the estimation of the 
channel impulse response by the channel impulse response 37. 



Figure 4 illustrates an improved version of the circuitry of Figure 1 . In Figure 4, the 
estimated symbol s is fed back to the channel impulse response unit 37 via a 
feedback decision generator 20 which generates a feedback decision D from the 
estimated symbol s. That feedback decision D is used to improve the estimate of 
channel impulse response generated by the CIR unit 37. An SIR block 29 receives 
the incoming signal samples x and the estimated signal u and calculates the signal 
to information ratio used to generate noise variance a 2 for the feedback decision 
generator 20. 

Figure 5 illustrates how the feedback decision D is used to improve the estimated 
channel impulse response. Figure 5 illustrates a CIR unit 37 which is capable of 
providing an estimated channel impulse response H using both fedback estimated 
data symbols and incoming pilot symbols. The circuit comprises first and second 
matched filters 10,12 which are each used to apply the relevant despreading code to 
the incoming digital signal x in . The reference code generator 22 generates three 
codes. A spreading long code LC is applied to both of the match filters 10,12. A 
spreading short code W p for pilot symbols is applied only to the first match filter 10 
and a spreading short code W d for data is applied only to the second match filter 12. 
The output signals 14,16 from the first and second match filters 10,12 respectively 
are applied to an averaging circuit 18. The averaging circuit receives the decision 
feedback D and applies that to the signals 14,16 from the matched filters to generate 
the estimated impulse response H. 



The averaging circuit 18 includes a memory 40 for holding previous channel impulse 
responses and a weighting/averaging unit 42 for taking an average of the earlier 
channel impulse responses and applying a weighting function according to the 
fedback decision value D. 

The correlation process which is used to establish the estimate of channel impulse 
response can be expressed as: 

H n (t)= I D n -C*(i)x(i+t)=Dn I C*{i)x(i+t) (Equation 1) 

i i 

In the above equation, i is the sampling time index along the received signal, t is the 
resulting sampling time index of the correlation process, x(i+t) is the received data 
sample for that index, C* is the conjugate of the spreading code and D n is the nth 
decision feedback where n is the time index of the impulse response. It will be 
apparent that the decision feedback D may be in respect of pilot symbols and/or data 
symbols. 

In use, an initial estimate for the channel impulse response H would be established 
by the CIR unit 37 on the basis only of known incoming pilot symbols in a particular 
time slot in a manner which is known perse. Alternatively, an initial estimate can be 
made using Equation (1) without the weighting function D. Once a preliminary 
estimate of the channel impulse response has been established, this is held in the 
memory 40 and supplied to the rake receiver 36 for subsequent signal processing of 
data symbols in that time slot. The first data symbol to be estimated is then subject 



to feedback to provide a feedback decision D so that a new estimate for the channel 
impulse response H can be determined using the procedure outlined above and as 
illustrated in Figure 5. The spreading code conjugate C* in equation 1 is derived 
from merging the conjugates of W d and LC (for data) or W p and LC (for pilot 
symbols). Then, an average is taken of the new estimate and the preliminary 
estimate. Thus, at each estimated symbol an improved channel impulse response is 
obtained for subsequent processing. The average can be taken coherently within 
any coherent period, that is part of a slot, one whole slot or a multiple of slots. 
Alternatively, a non-coherent average can be taken over any period. 

It will be apparent from Equation 1 that the correlation process can effectively be 
implemented without decision feedback as if all the transmitted data symbols are 
ones. The transmission feedback D acts as an average weighting factor to 
determine the estimated channel impulse response H. 

Another aspect of the implementation discussed herein is to determine the nature of 
the decision feedback D generated by the feedback decision generator 20. Figure 6 
illustrates the theory underlying one proposal which is to utilise a soft symbol 
decision feedback (SDFB) process by generating a soft feedback decision value D. 
Figure 6 illustrates a vector diagram of signal space. Vector A is a transmitted data 
signal. Vector A 1 is the received data. Vector A is the actual value we expect from 
the transmitted data signal since the rake receiver cannot be perfect. The received 
data symbol can be written as: 



A-a+n 
where n is the noise vector. 



(Equation 2) 



It is natural to assume that both the signal and the noise are zero mean. We have: 

E{a} =0 
E{n} =0 

E{A*}=E{a+n} =0 (Equation 3) 

In a binary (+1/-1) environment, we have also: 
E{AA*} =1 

E{A7V*} =E{(a+n)(a+n)*} 
=E{a a*}+ E{nn*} 

=E{a a*}+ a 2 (Equation 4) 



where a 2 is the noise variance. 



With these statistics of the signal, the soft output can be adjusted to satisfy the 
conditions for feedback if the noise variance is known beforehand. For a simple 
solution, we suppose that the signal loss is minimum and can be neglected: 

U-a| -> o 



Then equation (3) becomes: 



E{A'A'*} =E{a a*}+ a 2 
=E{AA*}+ a 2 



=1 + a 2 



(Equation 5) 



Using this analysis, the feedback decision generator can generate a soft value D 
satisfying the following requirements: 



As an alternative, a hard decision feedback may be utilised by performing a binary 
cut-off, that is to provide either +1 or -1 to the input of the channel impulse response 
unit for each estimated symbol S. 

The effectiveness of both solutions is illustrated in Figures 7 and 8. In figures 7 and 
8, the x axis represents signal energy per bit to noise ratio (E b /N 0 ) and the y axis 
represents the scale of acquisition probability - that is the likelihood that all signal 
propagation paths have been captured in the received signal as estimated using the 
C.I.R. Figure 7 is the simulation results with the best possible channel estimation, 
and Figure 8 is the simulation results in which the channel information of a slot is 
obtained only from the pilot symbols in the same slot. These simulations have been 



E{D} = 0 



(Equation 6) 



E{DD*}= 1+o 2 



(Equation 7) 



generated based on the MATLAB program included as Annexe 1. The simulation 
results compare the acquisition probabilities of the impulse response of different 
approaches of the decision feedback compared with no decision feedback. In both 
cases, it can be seen that the best results are given by a soft decision feedback for 
both pilot and data symbols. This is the top line in both simulations. The next line is 
for soft feedback decisions relating only to data symbols, that is without pilot symbol 
feedback being taken into account. 

The next line is for hard data feedback decisions in respect of data symbols only. 
The lowest line is where there is no decision feedback. 



ANNEXE 1 



MATLAB program 

% Soft Decision Feedback 

% ARIB specification 

% QZD/NTC - FEB. 98 - Oulu 

clear, rand( ' seed' , 0) , randn ( 1 seed 1 , 0) 

Dat_Weight = 0.5; SNR = -30:3:0; Fingers = 10; 

h= [0.7+j*0.1,0,0 f 0, -0.5- j*0. 3, 0,0, 0,0, 0.3+j*0. 5, 0,0, ... 

0.1-j*0.13,0,0, 0,0,0,0,0.3-j*0.1,0,0,0,0.5+j*0.5, ... 

0,0,0, -0.5- j*0. 3, 0,0, 0,0,0. 3 + j*0. 5, 0,0,0. l-j*0. 13, .. . 

0,0,0,0,0,0,0.3-j*0.1] ; 
h=h/sqrt(h*h') ; 

AcqError = zeros (size (SNR, 2) , 4) ; 

for kkk = 1: size (SNR, 2) 

snr = SNR (kkk) % SNR in dB 

for nnn=l:100 0 
DAT = (rand(l,40) >0.5) *2-l; 
% CODE 

CL_I= (rand (1,2560) >0.5)*2-l; CL_Q= (rand (1 , 2560 ) >0.5) *2-l; 
C = CL_I + j*CL_Q; 

w=l;for k=l:6;w=[w w;W -w] ; end; w0=w (1 , : ) ; wl=w (2 # : ) ; clear W; 
W0=[] ;W1=[] ; for k=l:40, W0= [W0 , wO] ; Wl= [Wl , wl] ; end 
CO = W0 .*C; CI = Wl .*C; 
% Spreading 

D=[]; for k=l:40, D- [D , DAT (k) *wl] ; end % Data 
P=ones(size(D) ) ; % PILOT 

xl = (D + j*P) .* C; 
% Channel 

x2 = conv(xl,h) /2; 
% Noise 

nA=l/10 A (snr/20) ;for k=l : size (x2 , 2) ; 

n(k) =nA*randn*exp(j*rand*2*pi) ;end 
x = x2 + n; 

%%%%%%%%%%%%%%%%%%% MF %%%%%%%%%%%%%%%%%%% 
IRPilot = zeros (39, 64); IRData = zeros (39, 64); 
for k=l:39, l=(k-l)*64; 
for m = 1:64 

IRPilot (k,m) = IRPilot (k,m) + x(l+m-l+ (1 : 64) ) *C0 (1+ (1 : 64) ) ' ; 
IRData (k,m) = IRData (k, m) + x(l+m-l+ (1:64) ) *C1 (1+ (1:64) )' ; 
end 

end 

IRPilot = -j*IRPilot/64; 
IRPilotHard=sum ( IRPilot ( 1 : 28 , : ) ) /28 ; 
IRData = IRData/64; 

%%%%%%%%%%%%%%%%%%% H< %%%%%%%%%%%%%%%%%%%%%%%%%% 
a = IRPilotHard . * conj (IRPilotHard) ; 
H = zeros (size (h) ) ; 
for k = 1: Fingers, findf = 0; 
while ( findf < 1) 

[n,m] =max(a) ; a (m) =0; 



if (m<=size(h,2) ) , H (m) -IRPilotHard (m) ; findf=2; end 

end 

end 

%%%%%%%%%%%%%%%%%%% decision %%%%%%%%%%%%%%%%%%% 
hl=find(abs(h)>0.01) ; 

h2=hl; % Idea 

%h2=f ind (abs (H) >0 . 01) ; % blind 

PilotSoft = zeros(l,28); DataSoft = zeros(l,40); 
for n=l: Fingers 

a = x(h2 (n) : size (CO , 2 ) +h2 (n) -1) .* con j (CO) * conj (H (h2 (n) ) ) ; 
for k=l:28 

PilotSoft (k) = PilotSoft (k) + imag (sum(a ( (k-1) *64+l :k*64) ) ) ; 
end 

a = x(h2(n) :size(Cl,2)+h2(n)-l) .* conj (CI) * conj (H <h2 (n) ) ) ; 
for k=l:40 

DataSoft (k) = DataSoft (k) + real (sum(a ( (k-1) *64+l :k*64) ) ) ; 
end 

end 

PilotSoft = PilotSoft/64; 
DataSoft = DataSoft/64; 
DataHard = sign (DataSoft ) ; 

%%%%%%%%%%%%%%%%%%% Averaging %%%%%%%%%%%%%%%%%%% 
a = 0; b = 0; c - 0; 

IRHard=zeros(l,64) ; IRSof t=zeros (1 , 64 ) ; IRSof tp=zeros (1, 64) ; 
IRSofter=zeros(l / 64) ; IRSof terp=zeros (1 , 64) ; 

for k=l:39 . 
IRHard=IRHard+IRData(k, : ) *DataHard (k) ; a=a+abs (DataHard (k) ) ; 
IRSof t=IRSof t+IRData (k, : ) *DataSof t (k) ; b=b+abs (DataSoft (k) ) ; 

end 

for k=l:28 

IRSof tp= IRSof tp+IRPi lot (k, : ) *PilotSof t (k) ; c=c+abs 
(PilotSoft (k) ) ? 

end , 
IRSoftp = IRSoft*Dat_Weight/b + IRSoftp * (l-Dat_Weight) /c; 
IRHard = IRHard*Dat_Weight/a + IRPilotHard* (l-Dat_Weight) ; 
IRSoft = IRSof t*Dat_Weight/b + IRPilotHard* (l-Dat_Weight) ; 
%************ acq probability ************ 
a = IRPilotHard .* conj (IRPilotHard) ; 
for k=l:size(hl,2) , [n,m] =max(a) ;h2 (k) =m ; a (m) =0 ;end 
c=0; for k=l:size(hl / 2) , f or n=l : size (hi , 2 ) 

if (h2(k)==hl(n)) , c=c+l; h2(k)=0; end 
end, end 
AErrNoDFB = size (hi , 2 ) -c ; 
a = IRHard . * conj (IRHard) ; 

for k-1: size (hi, 2) , [n,m] =max (a) ;h2 (k) =m ; a (m) =0 ;end 
c=0; for k=l:size(hl,2) ,for n=l : size (hi , 2) 

if (h2(k)==hl(n)), c=c+l; h2(k)=0; end 
end, end 
AErrHardDFB = size (hi, 2) -c; 
a = IRSoft .* conj (IRSoft) ; 

for k=l:size(hl,2) , [n,m] =max(a) ;h2 (k) =m ; a (m) =0 ; end 
c=0; for k=l:size(hl,2) , f or n=l : size (hi , 2 ) 



if (hi (k) ==hl (n) ) , c=c+l; h2(k)=0; end 
end , end 
AErrSof tDFB = size (hi , 2 ) -c ; 
a = IRSoftp . * conj (IRSoftp) ; 

for k=l : size (hi, 2) , [n,m] =max (a) ;h2 (k) =m;a (m) =0;end 
c=0; for k=l : size (hi, 2) , for n=l : size (hi, 2) 

if (h2 (k)==hl(n) ) , c=c+l; h2(k)=0; end 
end , end 
AErrSof tPDFB = size (hi, 2) -c; 
AcqError (kkk, : ) =AcqError (kkk, :)+,.. 
[AErrNoDFB , AErrHardDFB , AErrSof tDFB , AErrSof tPDFB] 
end % nnn 
end % kkk 



